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1.  Introduction 


An  effort  has  been  made  within  the  US  Anny  Research  Laboratory  (ARL)  to  extract 
quantitative  information  on  explosive  performance  from  high-speed  imaging  of 
explosions.  Explosive  fireball  surface  temperatures  are  measured  using  imaging 
pyrometry  (2 -color  2-camera  imaging  pyrometer;  full-color  single-camera  imaging 
pyrometer).  Framing  cameras  are  synchronized  with  pulsed  laser  illumination  to 
measure  fireball/shock  expansion  velocities,  enabling  calculation  of  peak  air-shock 
pressures.  Multicamera  filtering  at  different  wavelengths  enables  visualization  of 
light  emission  by  some  reactant  species  participating  in  energy  release  during  an 
explosion.  Measurement  of  incident  and  reflected  shock  velocities  is  used  to 
calculate  shock  energy  on  target. 

Results  of  these  measurements  are  used  to  construct  maps  of  temperature,  pressure, 
reactant  species,  and  shock  energy  on  target.  This  information  is  valuable  to 
evaluate  explosive  performance,  models  of  performance,  and  barriers  designed  to 
enhance  protection  and  survivability.  These  techniques  and  instruments  were 
developed,  in  part,  to  improve  productivity  by  lowering  testing  costs,  allowing  a 
single  event  to  yield  temperature,  pressure,  chemical  species,  and  perfonnance  data. 

2.  Pyrometry 


2.1  Background 

Pyrometry  is  the  method  of  estimating  temperature  of  incandescent  bodies  from 
standoff,  or  noncontact  methods.  Modem  optical  techniques  for  temperature 
measurement  began  in  the  early  1900s  using  the  disappearing-filament  technique, 
where  the  eye  of  an  operator  was  used  to  judge  when  the  red-filtered  incandescence 
from  a  sample  to  be  measured  matched  the  output  from  a  similarly  filtered, 
resistively  heated  filament  (the  filament  “disappeared”  when  the  2  images  were 
superimposed).1  By  1905,  the  techniques  had  advanced  to  the  point  of  publication 
of  a  National  Bureau  of  Standards  monograph  that  described  single  and  2-color 
optical  pyrometry  and  was  based  in  part  on  theoretical  work  on  gray  emitters  by 
Planck,  Wien,  and  others.2  The  commercial  application  of  optical  pyrometry  for  the 
next  50  years  relied  on  the  disappearing-filament  technique.  Accuracy  by  the  1960s 
approached  1  °C,  over  a  range  of  approximately  775-2,850  °C.3  With  the  advent  of 
solid-state  photodetectors  beginning  in  about  1960,  2-color4  or  “ratio”  optical 
pyrometry  began  replacing  disappearing-filament-based  instrumentation. 
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A  requirement  of  the  2-color  color  technique  is  that  variance  in  emissivity  with 
wavelength  is  small,  such  that  the  object  in  question  exhibits  gray  body 
behavior.5  Three-,  4-,  and  6-color  pyrometry  has  enabled  temperature  measurement 
of  materials  with  wavelength-dependent  emissivities,  provided  that  the  change  is 
continuous.6  In  2006,  Goroshin  et  al.  used  multiwavelength  pyrometers  to  measure 
temperatures  of  explosions  with  microsecond  time  resolution.7  Even  when 
accounting  for  changes  in  emissivity  with  wavelength,  the  spectral  region  to  be 
investigated  should  be  free  of  discrete,  gas  phase  emission.  Therefore,  each 
temperature  measurement  should  be  accompanied  by  a  full  emission  spectrum. 
Standoff  temperature  estimation  using  the  line  spectra  of  gas  phase  species  in 
emission  has  been  shown  to  be  effective  for  species  that  rapidly  reach  local  thennal 
equilibrium  with  their  surroundings  (e.g.,  some  atomic  emission).8  However,  fitting 
discrete  spectral  features  in  emission  can  be  risky  and  for  some  species  yields 
erroneous  values  of  temperature  (e.g.,  CO  emission  in  some  flames).9 

Extending  these  techniques  to  imaging  pyrometry  employing  light-sensitive  pixel 
arrays  began  in  the  1990s  with  work  on  thin  filament  pyrometry.10  With  the  advent 
of  fast  charged-coupled  device  (CCD)  cameras,  it  became  possible  to  image 
temperatures  and  to  follow  temperature  change  in  dynamic  events  (e.g.,  diesel 
engine  combustion).1113  Along  with  the  development  of  the  color  CCD  camera 
came  single-camera  pyrometry  that  employed  a  multiwavelength  technique  using 
the  R-G-B  Bayer-type  mask  for  pixel  filtering.14  Kuhn  et  al.  used  a  Nikon  D70  to 
temperature-map  an  atmospheric  pressure  flame.15  Densmore  et  al.  have  applied 
the  2-color  and  multiwavelength  imaging  techniques  to  explosions.16,17 

2.2  Theory 

The  principle  behind  most  pyrometric  techniques  is  an  application  of  Planck’s 
radiation  law  (Eq.  1),  where  U  is  light  intensity  at  wavelength  A,  h  is  Planck’s 
constant,  c  is  the  speed  of  light,  k  is  Boltzmann’s  constant,  si  is  the  emissivity  of 
the  material  (wavelength  dependent),  and  T  is  the  temperature.18 

Ii  =  2sihc2/A5  x  l/[exp(hc/AkT)  -  1]  .  (1) 

It  is  possible  to  use  Eq.  1  to  measure  the  temperature  of  a  body  with  one  detector. 
However,  this  is  impractical  as  the  intensity  of  light  falling  on  the  detector  is  highly 
optical-path  dependent  and  the  light  transfer  function  (how  the  optics  modify  light 
throughput)  must  be  carefully  characterized.19  To  mitigate  this  limitation,  the 
intensity  of  light  at  2  or  more  wavelengths  is  measured.  If  the  optical  paths  are 
similar,  the  transfer  function  can  be  removed  from  consideration.  The  ratio  of  light 
intensities  (Eq.  1)  at  2  wavelengths20  is  given  by  Eq.  2: 
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hi  / hi  =  [ sufa5  (exp(hc/fakT)  -  1)]  /  [ anfa5  (exp(hc/XikT)  -  1)]  .  (2) 

Because  the  emissivity  s  is  usually  not  known  for  the  body  under  study,  it  is  (often 
erroneously)  assumed  to  be  wavelength  independent  (gray  body  assumption)  and 
therefore  can  be  factored  out.21  As  the  equation  is  nonlinear,  2  approaches  have 
been  investigated  to  solve  Eq.  2  for  temperature.  The  first  is  to  assume  that,  at  the 
temperatures  of  interest  for  detonations,  the  exponential  tenn  is  large  compared  to 
1  (Wien’s  approximation22).  Dropping  one  from  the  numerator  and  denominator 
allows  for  an  analytical  solution.  Alternatively,  Eq.  2  is  kept  intact  and  a  numerical 
search  for  the  corresponding  temperature  is  performed.23 

2.3  Instrumentation 

In  work  described  here,  time-resolved  temperature  maps  of  detonations  of 
explosives  are  made  using  a  full-color  single-camera  pyrometer  where  wavelength 
resolution  is  achieved  using  the  Bayer-type  mask  covering  the  sensor  chip17  and  a 
2-color  imaging  pyrometer  employing  2  monochrome  cameras  filtered  at 
wavelengths  of  700  and  900  nm,  respectively,  (wavelength  regions  adjustable). 
Each  rig  operates  on  the  assumption  of  gray  body  behavior  described  previously, 
but  each  is  specific  to  the  type  of  explosive  being  investigated.  For  many  CHNO- 
based  explosives  (e.g.,  TNT  [C7H5N3O6],  the  fonnulation  C-4  [92%  RDX, 
CsEkNeOe]),  hot  detonation  products  are  mainly  soot  and  permanent  gases, 
presenting  an  approximation  of  a  gray  body  emitter.7,24  For  these  systems,  the 
single-camera  rig  may  be  appropriate.  For  metalized  explosives,  narrow-band  light 
emission  from  gas  phase  molecular  and  atomic  species  (e.g.,  AlO  near  484  nm,  BO2 
near  560  nm,  and  K  near  760  nm)  necessitates  the  use  of  the  2-camera  rig  to  make 
measurements  in  spectral  regions  free  of  discrete  features.  Additionally,  strong  C2 
or  CH  emission  from  nonsooting  explosive  fireballs  may  present  a  significant 
source  of  error.25  For  each  measurement  approach,  it  is  mandatory  to  measure  a 
time-resolved  emission  spectrum  during  the  event  to  ensure  the  absence  of  discrete 
emission  in  the  spectral  window  used  for  temperature  measurement.26 

For  each  pyrometer  rig  described  here,  framing  speeds  are  20,000-40,000  frames 
per  second  (fps)  at  a  resolution  of  approximately  400  x  500  pixels  with  an  exposure 
per  frame  of  one  to  tens  of  microseconds.  Each  system  is  temperature  calibrated 
using  a  standard  blackbody  source  (Omega  Engineering),  and  checked  for  accuracy 
using  an  air/acetylene  diffusion  flame.27 

A  schematic  of  the  ARL  imaging  pyrometer  rig  (showing  the  2-camera  rig  and  the 
full-color  rig  as  described)  is  shown  in  Fig.  1,  as  employed  for  temperature 
measurements  of  exploding  227-gram  (g)  spheres  of  the  explosive  formulation  C-4. 
Also  shown  in  this  figure  are  a  spectrograph  capable  of  measuring  time-resolved 
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emission  spectra  and  a  3 -color  spatially  integrating  pyrometer  used  as  a  check  of 
the  imaging  devices.  The  3-color  spatially  integrating  pyrometer  (700, 
850,  and  1,000  nm;  10-mn  bandpass)  can  provide  submicrosecond  time  resolution 
but  is  biased  toward  measuring  the  hottest  portion  of  an  emitting  medium  because 
of  the  T4  dependence  of  intensity.28 


Fig.  1  A  schematic  of  the  ARL  imaging  pyrometer  rig  as  employed  for  temperature 
measurements  of  exploding  spheres  of  the  explosive  formulation  C-4 

2.4  Results 


Figure  2  shows  a  sequence  of  temperature  maps  of  exploding  C-4  spheres  in  the 
near-field,  measured  using  the  single-camera  pyrometer.  These  images  were 
obtained  using  the  setup  shown  in  Fig  1.  Note  the  hot  outer  surface  of  the  fireball 
indicating  the  primary  shock  wave  location. 
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Fig.  2  Full-color  imaging  pyrometry  measurements  for  a  227-g  charge  of  C-4  at  63.5-cm 
standoff  from  a  metal-surface  table  used  to  simulate  an  armor  surface.  Time  sequence: 
a<b<c<d<e<f  with  a  time  interval  of  71  ps  between  frames.  X-  and  y-axes  labels  are  in 
millimeters.  Color  chart  bar  indicates  temperature  in  Kelvin. 

Figure  3  shows  a  frame  from  a  temperature  movie  of  a  fireball  produced  by 
detonation  of  a  1,000-g  TNT  charge  encased  in  a  boron/polymer  matrix.  Emission 
bands  from  BO2  near  546  mn  mandates  the  use  of  the  2-camera  rig.29  Regions  of 
the  image  in  Fig.  3  where  no  data  are  shown  were  either  saturated  in  one  or  both 
channels  or  had  intensity  levels  in  one  or  both  channels  below  the  camera  noise 
threshold. 
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Fig.  3  A  frame  from  a  temperature  movie  of  a  fireball  produced  by  a  1,000-g  TNT  charge 
encased  in  a  boron/polymer  matrix.  The  field  of  view  of  the  camera  is  approximately  2  m;  the 
frame  is  approximately  100  ps  after  initiation.  The  color  temperature  scale  is  in  Kelvin.  The 
field  of  view  is  approximately  1  m2.  Color  chart  bar  indicates  temperature  in  Kelvin. 

3.  Optically  Based  Pressure  Measurement 


3.1  Background 

Modern  optical  pressure  measurement  methods  often  rely  upon  a  light-sensitive 
material  placed  within  the  pressure  field.  Methods  include  pressure-sensitive 
fluorescence  quenching,30  pressure-induced  shifts  in  fluorescence  wavelengths,31 
pressure-induced  changes  in  interferometer  configurations,’2  pressure-sensitive 
changes  in  transmission  of  polarized  light  through  a  photoelastic  medium,33  fiber¬ 
optic-based  Mach-Zender  interferometers,34  and  tunable  diode  laser  measurement 
of  pressure-induced  changes  in  absorption  line  profiles  of  gases.35  Standoff  optical 
pressure  measurements  relating  shock  velocity  to  peak  shock  pressure  can  trace 
their  origins  to  the  experiments  of  Hooke,  Huygens,  Mach,  and  others,  studying 
heat  and  flow  visualization  via  changes  in  refractive  index  of  air,36  and 
simultaneous  development  of  high-speed  photography.  Shock-wave  imaging  using 
a  short  pulse  light  source  in  a  schlieren  system  was  first  demonstrated  by  Mach  and 
Salcher37  in  1887.  In  the  early  twentieth  century,  the  schlieren  technique  developed 
by  Mach  became  the  most  important  tool  for  visualizing  flow  during  ballistic 
events.  In  1929,  the  first  high-speed  schlieren  movies  (1  million  fps  for  24  frames) 
were  demonstrated  by  Cranz  and  Schardin.38  The  development  of  high-speed 
photography  by  the  1940s  focused  on  high-rotation-speed  film39  or  rotating  mirrors 
or  prisms.40  Schlieren  imaging  became  the  most  important  tool  for  flow 
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visualization  through  the  1970s,  mainly  through  work  related  to  supersonic  flight 
and  rocket  design.36  Remote  measurements  of  peak  pressures  in  explosions  usually 
rely  on  measuring  shock  wave  velocities  and  use  Rankine-Hugoniot  relations  to 
convert  shock  Mach  numbers  to  pressures.41  Improvements  in  optics  and  recording 
techniques,  led  by  the  work  of  Settles  et  ah, 36  enable  highly  sensitive  measurements 
of  shock  position  following  explosions  and  of  changes  in  refractive  index  behind 
the  shock.42 

Tracking  explosively  generated  shock-wave  position  using  lasers  was  first  reported 
in  the  late  1980s.43  Pulsed  copper  (Cu)-vapor  lasers,  synchronized  to  framing 
cameras  at  repetition  rates  up  to  10  kHz,  were  first  reported  in  1987. 44  In  2002, 
Frost  et  ah,  using  a  combination  of  high-speed  imaging  and  piezo-resistive  gauges, 
reported  acceleration  of  the  shock  front  by  metal  particle  combustion,  one  of  the 
first  investigations  of  the  explosive  near-field.45  In  2005,  a  20-kHz  Cu  vapor  laser 
(5 10-nm  output)  was  synchronized  to  a  high-speed  filtered  camera  to  obtain  fireball 
expansion  and  shock  velocities.  Peak  pressures  calculated  from  movies  of  the  shock 
wave  agreed  reasonably  well  with  pressures  reported  by  mechanical  transducers.46 

3.2  Theory  and  Instrumentation 

The  pressure  field  of  a  detonation  event  is  often  measured  at  discrete  points  in  space 
where  gauges  are  located.  Many  piezoelectric-based  gauge  technologies  suffer 
from  both  thermal  effects  and  fragility.47  49  This  often  limits  conventional  pressure 
measurements  to  larger  standoff  distances  (j£5  charge  diameters).  The  optical 
technique  used  here  overcomes  these  limitations,  but  is  limited  in  that  it  cannot 
provide  impulse  or  pressures  within  the  fireball.  However,  data  from  the  single  or 
2-camera  pyrometer  can  capture  the  fireball  and  subsequent  air  shock  structure  of 
the  detonation  event,  allowing  simultaneous  temperature  and  pressure 
measurement.  Utilizing  the  Rankine-Hugoniot29  relationship  between  shock 
velocity  and  shock  pressure  for  a  spherical  shock  expanding  into  an  ideal  gas,  where 
Ps  is  the  shock  pressure,  Pa  is  ambient  pressure,  and  M  is  Mach  number  based  on 
measured  shock  velocity,  the  information  from  the  pyrometer  can  be  used  to  map 
the  pressure  field  of  the  explosion  (Eq.  3). 

Ps/Pa  =  7[M2-l]/6.  (3) 

To  calculate  shock  pressures,  it  is  assumed  that  at  early  times  (time  after  detonation 
of  1  ms  or  less)  the  shock  wave  is  coincident  with  the  outer  edge  of  the  expanding 
fireball.  At  later  times,  when  the  air  shock  has  separated  from  the  fireball  contact 
surface,  the  shock  wave  position  is  used  to  calculate  shock  velocities.50  Figure  4 
shows  an  image  sequence  of  center-detonated  300-g  exploding  C-4  spheres 
measured  using  a  high-speed  framing  camera  (Cordin  Model  570).  Exposure  time 
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per  picture  is  300  ns,  and  the  time  between  each  exposure  is  approximately  1  ps. 
For  longer  times  and  slower  expansion  velocities,  a  second  digital  framing  camera 
(Photron  SA-5)  was  used  to  record  the  fireball  expansion.  Consecutive  images  are 
used  to  determine  shock  front  position  versus  time,  and  a  peak  pressure  calculated 
from  measured  shock  velocities,  employing  the  Rankine-Hugoniot  relations 
(Eq.  3). 


Fig.  4  An  image  sequence  of  center-detonated  exploding,  300-g  C-4  spheres  measured 
using  a  high-speed  framing  camera.  Exposure  time  per  picture  is  300  ns,  and  the  time  between 
each  exposure  is  approximately  1  ps. 

3.3  Results 

Figure  5  shows  the  result  of  peak  shock  pressures  calculated  using  data  from  the 
high-speed  framing  cameras  (Cordin  570  and  Photron  SA-5)  and  results  calculated 
using  the  blast  peak  pressure  simulator  CONWEP.51  Agreement  is  reasonable 
overall.  For  this  multiple  camera  rig,  error  in  the  measurements  is  believed  to  be 
less  than  10%  at  all  distances  from  charge  center,  with  the  main  sources  of  error 
being  the  duration  of  exposure  (need  for  even  faster  cameras  at  very  early  times) 
and  deviation  from  ideal  gas  behavior  at  higher  Mach  number.29 
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Fig.  5  Optically  based  peak  shock  pressure  measurements  using  the  framing  camera 
(Cordin  Model  570)  and  the  streak  camera  rig,  and  predictions  based  upon  CONWEP.  The 
data  is  from  the  images  shown  in  Fig.  4. 

4.  Chemical  Species  Imaging 


4.1  Background 


Several  metals  and  metalloids  (e.g.,  A1  [aluminum],  B  [boron],  and  Si  [silicon]) 
may  be  considered  for  addition  to  propellants  and  explosives  to  increase  energy.40 
For  example,  Eqs.  4-7  show  the  combustion  energies  (enthalpy  of  combustion, 
AHC),  in  kilojoules  (kJ)  per  mole  (mol)  and  kJ/g  of  fuel  for  the  oxidation  reactions 
of  B,  Al,  and  Si,  respectively.  The  combustion  enthalpy  of  trinitrotoluene  (TNT, 
C7H6N3O6,  Eq.  4)  is  also  shown  for  comparison.52 

4B  +  302  — >  2B2O3  (i)  AHC  =  -626.8  kJ/mol  AHC  = -57.97  kJ/g  (4) 

4A1  +  302  — >■  2Ab03  (i)  AHC  =  -810  kJ/mol  AHC  = -30.03  kJ/g  (5) 

Si  +  O2  — >  SiC>2  (vitreous)  AHC  = -910.7  kJ/mol  AHC  =  -32.4  kJ/g  (6) 

2  C7H6N3O6  +  1 IO2  — *  14C02  +  6H2O  +  3N2 


AHC  =  -34 1 0  kJ/mol  AHC  =  -15.1  kJ/g  (7) 


Because  of  the  significant  gain  in  combustion  energy  (weight  basis,  relative  to 
TNT)  predicted  by  combustion  calculations,  B  has  received  considerable  attention 
as  an  additive  to  explosive  and  propellant  mixtures.53  However,  practical 
thermodynamic  and  physical  properties  inhibit  realization  of  this  ideal.54  In  a 
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review  article,  Yeh  and  Kuo53  summarized  2-stage  boron-particle  combustion  and 
measured  burn  times.  Stepwise  reactions  occurring  in  B  oxidation  include  55 

B  +  O  =  BO  (8) 

BO  +  O  =  B02  (9) 

Reactions  shown  in  Eqs.  8  and  9  are  of  special  interest  for  studies  of  B  oxidation56 
and  the  work  described  here  because  excited  states  of  the  product  species  are 
responsible  for  the  strong  visible  (green)  light  emission,  and  offer  the  possibility  of 
studying  their  emission  as  a  function  of  ignition  stimuli.  Johns  has  shown29  that  the 
band  structure  seen  in  emission  and  fluorescence  during  B  combustion  arises  from 
sequences  of  2  electronic  transitions: 


A2nu  -  x2ng 

(10) 

b2e+u  -  x2ng 

(ii) 

A  goal  of  this  work  is  to  map  the  onset  of  B  oxidation  during  an  explosive  event. 

4.2  Theory  and  Approach 

Figure  6  shows  the  wavelength  spectrum  in  emission  of  a  burning  1 : 1  (by  weight) 
mixture  of  a  military-specification-compliant  (MIL-STD-961)57  B/KNO3  igniter 
formulation.  Arrows  show  locations  of  strongest  observed  emission  bands  of  the 
BO2  molecule.  Off-scale  features  near  589  and  760  mn  are  emission  from  Na 
(sodium)  and  K  (potassium)  atoms,  respectively.  To  obtain  this  spectrum, 
approximately  250  mg  of  the  igniter  fonnulation  was  placed  on  a  resistively  heated 
NiCr  (nickel  chromium)  wire  ribbon.  When  ignited,  visible  light  emission  was 
measured  using  an  Ocean  Optics  HR2000+  spectrograph  for  an  integration  time  of 
100  ms  after  trigger.  The  spectrograph  was  externally  triggered  by  a  Si  photodiode 
(New  Focus)  coupled  to  a  Stanford  Research  Systems  DG535  Pulse/Delay 
generator. 
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Fig.  6  Emission  from  a  burning  1:1  (by  weight)  powder  mixture  of  a  milspec  B/KNO3 
(MIL-STD-9617)  igniter  formulation.  Spectral  lines  from  gaseous  BO2  emission  (arrows) 
appear  superimposed  upon  particle  incandescence.  Not  corrected  for  intensity  response  of  the 
spectrograph. 

4.3  Instrumentation 

The  emission  spectrum  of  burning  B/KNO3  shown  in  Fig.  6  was  used  as  a  guide  to 
devise  a  scheme  for  imaging  emission  from  BO2  following  ballistic  initiation. 
Early-time  imaging  of  metal  oxidation  products  was  suggested  by  Dreizin  for  AlO 
emission  from  aluminized  explosives,  and  modified  here  for  imaging  of  BO2 
emission  (2012  December  personal  communication  with  EL  Dreizin, 
unreferenced).  Essentially,  2  light  sensors  (cameras),  each  filtered  over  a  narrow 
wavelength  region,  observe  an  event  over  the  same  line  of  sight.  The  wavelength 
regions  are  selected  so  that  one  region  coincides  with  a  known  emission  band  of 
the  molecule  of  interest.  The  second  region  is  selected  so  that  it  sees  only  particle 
incandescence  (gray  body  emission).36  The  light  signals  are  then  ratioed  to 
minimize  the  intensity  from  the  background  incandescence  (subtraction  gave  a 
qualitatively  similar  result).  For  imaging  BO2  emission,  the  light  sensors  were  2 
Phantom  V7.3  monochrome  cameras  (Vision  Research,  Inc.)  looking  through  a 
glass  plate  beam  splitter  (Edward  Scientific)  along  a  common  line  of  sight.  The 
cameras  were  run  from  the  same  time  base  (slave  and  master  setup)  using  identical 
lenses,  apertures,  frame  rate  (typically  40  kHz),  and  exposure  duration  (typically 
1  per  frame  rate).  For  measurements  of  BO2  emission,  one  50-mm-diameter  filter 
had  a  bandpass  (10  nm)  centered  at  546  nm,  and  the  other  50-mm-diameter  filter 
had  a  bandpass  (10  nm)  centered  at  700  nm.  Selection  of  these  optical  regions  was 
based  upon  the  spectrum  shown  in  Fig.  6.  Filters  were  obtained  from  Aurora 
Optical,  Inc.  The  sensitivity  of  the  cameras  at  wavelengths  of  546  and  700  nm  are 
approximately  equal36;  however,  no  correction  was  made  to  account  for  absolute 
differences  in  sensitivity. 
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4.3.1  2-Camera  Image  Registration 

Although  an  effort  was  made  to  adjust  the  cameras  to  a  common  line  of  sight,  it 
was  found  that  a  perfectly  matched,  hardware-based  pixel-by-pixel  overlap  from 
the  2  cameras  was  never  actually  attainable.  A  software  program  was  written  in  the 
computer  program  MATLAB  to  obtain  a  calibration  image,  create  a  transfer  matrix, 
and  then  apply  this  matrix  pixel-by-pixel  to  the  slave-camera  images.  Using  this 
software  program,  a  fresh  camera  placement,  and  fully  aligned,  registered  system 
could  be  achieved  in  approximately  10  min. 

4.3.2  Particle  Temperature  Measurements 

A  400-pm-diameter  Si-core-Si-clad  optical  fiber  (Aurora  Optics,  Hanover,  NH) 
was  inserted  into  the  top  of  the  test  chamber  and  the  distal  end  coupled  to  an  Ocean 
Optics  HR4000  visible  wavelength  spectrograph.  As  a  check,  light  launched  into 
this  optical  fiber  at  the  spectrograph  end  illuminated  approximately  the  full  volume 
imaged  by  the  2-camera  rig.  Therefore,  any  light  emitted  during  the  ballistic  event 
within  the  chamber  would  be  sampled  by  the  spectrograph.  Wavelength-resolved 
spectra,  triggered  by  the  gas  gun  firing,  were  recorded  at  15 -ms  intervals  during  the 
event  (timing  limited  by  the  spectrograph).  To  measure  temperature,  the  Ocean 
Optics  HR-4000  spectrograph  was  calibrated  for  intensity  response  using  an  Omega 
Inc.  Model  BB-4A  calibrated  blackbody.  Corrected  intensities  reported  at  850  and 
990  nm  were  used  to  calculate  a  gray-body  temperature  according  to  the  Planck 
function.28  This  wavelength  region  was  selected  because  it  is  free  from  discrete 
spectral  emission,  some  spectra  were  saturated  from  650  to  800  nm,  and  it  was 
decided  to  use  a  common  wavelength  range  for  all  temperature  calculations.  A 
check  of  the  temperature  measurement  technique  using  emission  from  an 
acetylene/air  diffusion  flame  gave  reasonable  results  (1,800  K  outer  soot 
temperature).27 

4.3.3  Light-Gas  Gun 

The  light-gas-gun  facility  at  General  Sciences  Inc.  (GSI,  Souderton,  PA)  was  used 
to  investigate  the  impact-velocity  dependence  of  ignition  of  the  B/KNCfi-filled 
polyethylene  spheres.  The  gas  gun  consists  of  a  high-pressure  helium  tank,  a  dump 
valve  that  exhausts  the  high-pressure  gas  to  the  gun  breech,  and  a  25-mm-diameter, 
unrifled,  4.9-m-long  gun  barrel.  The  projectile,  equipped  with  a  sabot,58  can  achieve 
velocities  approaching  1,200  m/s.  Precise  velocities  are  measured  photonically 
following  stripping  of  the  sabot  as  the  projectile  exits  the  barrel. 
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KNC>3  powder  (50- urn  particle  diameter)  and  B  powder  (30-uin  particle  diameter) 
were  obtained  from  Fisher  Scientific.  The  hollow  25-mm-diameter  polyethylene 
spheres  were  manufactured  in-house  at  GSI. 

Following  velocity  measurement,  the  spherical  projectile  travels  0.5  m  through  free 
space  and  then  enters  a  50-mm-diameter,  burst-disk-equipped  (0.18-mm  thickness, 
A1  sheet)  aperture  at  the  near  end  of  a  2-m3  catch  tank.  The  distal  end  of  the 
cylindrical  catch  tank  is  equipped  with  30-mm-thick  Pyrex  windows  that  allow 
visualization  of  the  projectile  impact  with  a  steel  anvil  at  the  far  interior  wall  of  the 
tank.  A  schematic  of  the  gas-gun  catch-tank  assembly  is  shown  in  Fig.  7,  along  with 
a  color  camera  image  (camera  mounted  to  an  optical  port  on  top  of  the  chamber)  of 
an  impact  event  using  the  B/KNCb-filled  polyethylene  projectile. 


Burst 
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Fig.  7  Schematic  of  the  gas-gun  catch-tank  assembly.  The  photograph  of  the  impact  event 
is  a  single  frame  from  a  color  camera  mounted  at  the  top  of  the  catch  tank  assembly. 

The  field  of  view  of  the  2-camera  rig  was  approximately  15  cm.  The  spherical 
polyethylene  projectile  weight  was  1 .25  g,  the  weight  of  the  1:1  by  weight  B/KNO3 
fill  was  1  g  (±  0.05  g),  and  total  projectile  weight  per  shot  was  approximately  2.25  g. 
The  2-camera  imaging  rig  was  positioned  to  enable  side-on  imaging  of  the 
B/KNCfi-filled  polyethylene  sphere  impacting  the  steel  anvil  at  the  rear  interior  of 
the  catch  chamber.  The  focal  plane  of  the  2-camera  rig  contained  the  projectile  shot- 
axis.  Figure  8  is  a  photograph  of  the  catch  chamber  with  the  2-camera  rig  in  position 
at  the  experimental  facility  at  GSI. 
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Fig.  8  The  2-camera  rig  in  position  at  the  side  port  of  the  catch  tank  prior  to  a  gun  firing. 
Note  the  color  camera  at  the  top  port  of  the  catch  tank.  The  gas  gun  is  positioned  to  the  left  of 
the  chamber  out  of  the  field  of  view  of  the  photo.  The  white  plastic  fixtures  in  front  of  each 
camera  lens  contain  the  optical  filters. 

4.4  Results 

The  monochrome  cameras  in  the  2-camera  rig  were  filtered  at  546  and  700  mn 
(each  using  a  10-nm  bandpass  width).  An  examination  of  the  emission  spectrum  of 
burning  B/KNO3  (Fig.  6)  and  comparison  with  measured  spectral  emission 
wavelengths  for  combusting  B/KNO3  mixtures59  shows  that  the  emission  near 
546  mn  is  a  combination  of  incandescing  hot  particles  and  gas  phase  BO2  emission, 
and  the  emission  near  700  mn  is  likely  due  only  to  incandescing  hot  particles. 

After  recording  a  movie  from  each  synchronized  camera,  and  registering  the  images 
as  described  previously,  data  were  prepared  by  doing  a  pixel-by-pixel  ratio  of 
intensities  (I546nm/l700nm)  of  each  corresponding  frame  from  each  camera.  If  the  ratio 
result  was  less  than  1,  indicating  higher  intensity  near  700  mn,  the  pixel  was  given 
a  red  color.  If  the  ratio  result  was  greater  than  1 ,  indicating  a  higher  intensity  from 
BO2  emission,  the  pixel  was  given  a  green  color.  Each  frame -to-frame  ratio 
produced  a  frame  of  a  new  movie  (Fig.  9),  consisting  of  green  and  red  pixels  only. 
The  goal  was  to  obtain  a  simple  measure  of  how  the  impact  velocity  influenced  the 
ratio  of  excited  BO2  molecules  to  incandescent  particles. 
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Fig.  9  A  series  of  processed  images,  time  after  impact,  that  are  part  of  an  emission  movie 
to  image  B  oxidation/BCh  emission.  The  projectile  enters  at  the  left  of  the  figure;  the  anvil  is 
at  the  right  of  the  figure.  Red  indicates  thermal  emission  (incandescence)  near  a  wavelength 
of  700  nm  has  higher  overall  intensity;  green  indicates  emission  from  incandescence  plus  BO2 
emission  near  a  wavelength  of  546  nm  has  higher  overall  intensity.  The  x  and  y  scales  are  pixel 
indices.  Each  pixel  corresponds  to  approximately  100  pm  in  length  at  the  focal  plane  of  the 
2-camera  rig. 

To  evaluate  whether  impact  velocity  has  a  real  influence  on  time  to  appearance  of 
BO2  emission  using  an  arbitrary  fixed  ratio  (I546nm/l700nm  =  1),  it  is  necessary  to 
measure  the  temperature  of  the  hot  particles.  Figure  10  shows  a  plot  of  average 
particle  (incandescence)  temperature  (850  nm,  990  nm)  during  time  of  peak  BO2 
emission,  and  time  to  onset  of  BO2  emission  ((I546nm/l700nm  >1)  versus  impact 
velocity  for  impact  experiments  of  B/KNC>3-containing  polyethylene  spheres.  Shots 
with  polyethylene  spheres  that  did  not  contain  the  B/KNO3  mixture  did  not  exhibit 
any  of  the  discrete  emission  here  attributed  to  B  oxidation/BCF  emission. 
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Fig.  10  A  plot  of  particle  incandescent  temperature  and  time  to  onset  of  gaseous  BO2 
emission  (i.e.,  ignition  delay)  vs.  impact  velocity  for  B/KNCb-containing  polyethylene  spheres 

4.5  Discussion 

Figure  10  shows  that  the  average  measured  particle  temperature  for  all  experiments, 
during  the  time  of  most  intense  BO2  emission,  is  near  1,400  K  (±100  K).  At  this 
temperature  the  intensity  of  BO2  gas  emission/particle  incandescence  near  546  nm 
exceeds  that  of  particulate  emission  intensity  alone  near  700  nm.  In  their  review, 
Yeh  and  Kuo  estimate  the  ignition  of  B  in  air  to  occur  at  approximately  1,900  K 
and  in  oxygen  at  approximately  1,580  K.53  No  data  are  given  for  ignition  of 
B/KNO3  mixtures.  The  data  shown  in  Fig.  10  indicate  that  onset  of  BO2  emission 
near  546  nm  in  the  experiment  described  here  is  a  direct  function  of  particle 
temperature  and  has  a  more  complicated  dependence  upon  impact  velocity.  That  is, 
once  the  powder  mixture  of  B  and  KNO3  reaches  approximately  1,400  K,  BO2, 
emission  becomes  intense. 

The  data  shown  in  Fig.  10  suggest  that  the  imaging  method  described  here  may  be 
of  use  in  investigating  when  a  metal  additive  participates  in  energy  release  in  an 
energetic  material  formulation.  The  impact  velocity  versus  onset  of  reaction  time 
data  shown  in  Fig.  10  suggest  that  at  low  velocity,  there  is  an  ignition  delay  during 
which  the  particle  is  heated  to  approximately  1,400  K  by  oxidation  reactions,  and 
the  ignition  delay  decreases  as  the  projectile  velocity  increase  to  800  m/s.  At  higher 
velocities,  there  is  an  apparent  rise  in  ignition  delay  with  impact  velocity.  It  is 
conceivable  that  the  apparent  increase  in  ignition  delay  time  at  impact  velocities 
above  around  850  m/s  may  be  caused  by  the  hottest  particles  leaving  the  camera 
field  of  view  as  these  particles  are  expected  to  have  high  velocities  relative  to 
particles  produced  at  lower  impact  velocities.  This  suggests  a  sweet  spot,  or  an 
optimum  impact  velocity  for  minimizing  time  to  onset  of  B  oxidation  with  these 
experimental  parameters  for  the  field  of  view  achieved  with  the  2-camera  rig. 
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Additionally,  although  it  is  implied  that  the  onset  of  green  pixels  in  the  individual 
images  indicates  ignition  delay  times,  the  onset  of  green  pixels  actually  means  that 
the  intensity  of  emission  at  546  nm,  at  this  time  after  impact,  exceeds  that  of 
incandescent  emission  at  700  nm.  For  the  data  shown  in  Fig.  10,  the  projectile 
velocity  at  impact  is  believed  to  have  an  error  approaching  10%,  the  uncertainty  in 
time  to  BO2  emission  is  believed  to  be  approximately  100  ps,  and  the  error  in 
calculated  temperature  is  believed  to  be  approximately  100  K,  based  upon 
calibration  measurements  using  emission  from  an  acetylene/air  flame. 

5.  Blast  Energy  on  Target 


5.1  Introduction 

High-speed  wavelength-filtered  imaging,  in  combination  with  pulsed  laser 
illumination,  has  been  used  previously45  to  image  shock  waves  near  the  expanding 
detonation  products.  Recently,  a  new  effort  was  initiated  to  measure  blast  energy 
on  target.  This  program  measures  incoming  and  reflected  blast  waves,  determines 
the  energy  in  each,  and  from  the  difference  determines  blast  energy  imparted  to  a 
target.  To  accomplish  this,  an  enhancement  in  fidelity  of  blast  wave  imaging  was 
needed. 

5.2  Instrumentation 

Figure  1 1  shows  a  schematic  of  a  low-parallax  high-brightness  imaging  rig  (LP- 
HBI)  developed  for  this  effort,  which  is  used  to  image  incoming  and  reflected  shock 
waves.  A  method  to  measure  the  blast  energy  deposited  on  the  target  is  to  calculate 
the  initial  energy  per  gram  of  shocked  air  in  the  incident  shock  wave  and  the  energy 
per  gram  in  the  reflected  shock  wave  immediately  before  and  after  the  encounter 
with  the  visor.  These  energies  are  calculated  using  only  the  peak  pressure  of  the 
respective  waves.  The  difference  is  assumed  to  be  the  energy  absorbed  by  the  target 
at  the  peak  incident  pressure  and  is  expressed  as  a  percentage  of  the  peak  incident 
energy  per  gram.  This  percentage  is  assumed  to  apply  to  the  total  impulse  of  the 
incident  wave. 
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~15  m 


Fig.  11  The  low-parallax  high-brightness  imaging  rig.  Low  parallax  is  achieved  by  placing 
the  laser  directional  mirror  at  the  center  of  the  camera  lens. 


5.3  Imaging 

Figure  12  shows  a  4-image  sequence  of  a  shock  wave  being  reflected  by  a  head 
mannequin  incorporating  personnel  protective  equipment.  For  the  discussion  that 
follows,  a  subscript  of  0  indicates  a  property  of  ambient  air,  a  subscript  of  1  or  2 
indicates  a  property  of  the  incident  or  reflected  shock,  respectively.  The  images 
shown  in  Fig.  12  were  recorded  at  16,000  fps,  for  an  interframe  time  of 
approximately  67  us,  at  an  exposure  time  of  1  ps.  Laser  pulse  duration 
(synchronized  to  the  camera)  was  10  ns  per  pulse  with  a  laser  energy  per  pulse  of 
approximately  1  mJ.  A  600-g  pentolite  (1:1  mixture  by  weight  of  TNT  and  PETN 
(C5FI8N4O12)60  sphere,  center  of  mass  2  m  from  the  visor  midpoint,  was  used  to 
provide  the  blast. 
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Fig.  12  A  shock  wave  being  reflected  by  a  mannequin.  Peak  pressure  at  the  visor  was 
approximately  69.5  psi.  Charge  was  660-g  pentolite,  2  m  from  center  of  face  shield. 

5.4  Results 

To  calculate  the  energy  imparted  to  the  target,  the  required  inputs52  are  as  follows: 

•  Measured  incident  (Ui)  and  reflected  (U2)  shock  velocities 

•  Hugoniot  for  air 

The  incident  and  reflected  shock  velocities  were  calculated  from  sequential  Photron 
SA-5  images  (Fig.  12  and  Table  1).  Time  of  images  in  Fig.  12  is  a<b<c<d.  The  air 
Hugoniot61  can  be  expressed  as  a  relation  between  the  shock  (Ui)  and  particle  (ui) 
velocities  in  meters  per  second  for  air  at  a  density  (po)  of  1 .29  kg/m3  or  specific 
volume,  Vo,  of  0.775  m3/kg. 

Ui=  0.2425  +  1.1231ui.  (12) 
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The  incident  shock  overpressure  may  be  given  by52 

Pi  -  Po  =  poUiui  =  (1.29  kg/m3)(705.2  m/s)(41 1.9  m/s)  =  374708.7  Pa  .  (13) 

Units  of  pressure  are  Pascals  (Pa),  Po=  10 1,325  Pa  and  Pi  =  absolute  shock  pressure. 
Using  these  relationships,  the  incident  peak  pressure  on  the  visor  in  Figure  12b  is 
approximately  476,033.7  Pa,  or  68.6  psi. 

The  initial  specific  volume,  Vo=  1/(0.00129  kg/m3)  =  0.775  m3/kg.  The  shocked 
specific  volume,  Vi,  may  be  given  by62 

Vi=  0.775(Ui-ui)/Ui  =  0.322  m3/kg  .  (14) 

From  Cooper,52  the  energy  equation  for  a  shock  moving  into  unshocked  air  is 

E2  =  l/2(Pi+Po)(Vo-Vi)  +  l/2(ui)2  .  (15) 

For  the  incident  shock  shown  in  Fig.  12b,  this  yields  a  specific  volume  shock  energy 
of  2 1 5 .6  kiloj  Joules  per  kilogram  (kJ/kg).  A  summary  of  physical  parameters  of  the 
shock  prior  to  the  visor  is  shown  in  Table  1 . 


Table  1  Physical  parameters  of  the  shock  prior  to  visor 


Ui 

Ul 

Vo 

Vi 

Po 

Pi 

Ei 

(m/s) 

(m/s) 

(m3/kg) 

(m3/kg) 

(Pa) 

(Pa) 

(kJ/kg) 

705.2 

411.9 

0.775 

0.322 

101,325 

476,033.7 

215.6 

After  encountering  the  visor,  the  portion  of  the  incident  shock  that  is  reflected  will 
travel  back  through  preshocked  air.  Relative  to  ambient  conditions,  this  air  has  a 
higher  density  and  a  particle  velocity,  ui,  caused  by  the  incident  shock.  The  particle 
velocity  in  the  preshocked  air  (ui)  is  added  to  the  observed  reflected  shock  velocity 
(Ur)  to  get  the  true  propagation  velocity  of  the  reflected  wave  at  the  moment  of 
reflection  into  the  shocked  air,  U2. 

Ur  +  ui=U2  =  366.7  m/s  +  41 1.9  m/s  =  778.6  m/s  .  (16) 

From  the  full  image  sequence  of  the  experiments,  it  was  ascertained  that  there  was 
no  perceptible  motion  of  the  visor  after  impact,  indicating  that  the  particle  velocity 
of  the  air  was  essentially  zero  behind  the  reflected  shock.  The  change  in  particle 
velocity,  dui,  is  then  equal  to  ui. 
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The  reflected  shock  travelled  into  air  with  a  specific  volume  imparted  by  the 
incident  shock.  The  specific  volume  of  air  behind  the  incident  shock,  V2,  is  given 
by 

V2  =  (Vi)(Ur/U2)  =  (0.322  m3/kg)(366.7  m/s  /  778.6  m/s)  =  0.152  m3/kg.  (17) 

The  reflected  shock  overpressure  is  then  calculated  as  in  Eq.  13. 

P2-P1  =  P1U2U1  =  (l/Vi)(U2)(dui)  =  (3.11  kg/m3)(778.6  m/s)(41 1.9  m/s) 

=  997,393.6  Pa.  (18) 

The  reflected  energy  per  gram,  E2,  is  equal  to  the  change  in  internal  energy  per 
gram  plus  kinetic  energy  per  gram. 

E2  =  1/2(P2+Pi)(V  1-V2X 1 06)  +  l/2(ui(105))2  (kJ/kg).  (19) 

E2  =  1/2(1949461  Pa)(0.17  m3/kg)  =165.7  kJ/kg.  (20) 

A  summary  of  physical  parameters  of  the  shock  after  reflection  by  the  visor  is 
shown  in  Table  2.  This  analysis,  although  preliminary,  indicates  that  for  the  target 
assembly  shown  in  Fig.  12  (mannequin  plus  visor),  the  reflected  shock  contained 
77%  of  the  energy  of  the  incident  shock.  Therefore,  the  target  assembly  experienced 
approximately  23%  of  the  initial  shock  energy.  Figure  13  shows  a  plot  of  incoming 
versus  reflected  shock  velocity  for  the  shock/mannequin  interaction  shown  in 
Fig.  12. 


Table  2  Physical  parameters  of  the  shock  after  reflection  by  visor 


Ur 

du 

Ul 

U2 

Vi 

V2 

Pi 

P2 

e2 

(m/s) 

(m/s) 

(m/)c 

(m/s) 

(m3/kg) 

(m3/kg) 

(Pa) 

(Pa) 

(kJ/g) 

366.7 

411.9 

411.9 

0 

0.322 

0.152 

476,033.7 

1473,427.3 

165.7 
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Fig.  13  Incoming  vs.  reflected  shock  velocity  for  Fig.  8. 

6.  Summary  and  Conclusions 

The  work  described  here  indicates  that  high-speed  high-fidelity  imaging 
incorporating  optical  filtering  and  signal  processing  techniques  can  make  a 
significant  contribution  to  analysis  of  the  behavior  of  high  explosives.  The  data 
generated  using  these  techniques  have  applicability  to  explosives  and  munitions 
formulation  and  design  and  to  improving  soldier  protection.  The  ability  to 
simultaneously  create  temperature,  peak  pressure,  and  chemical  species  maps  of 
explosions  offers  experimental  measurements  that  rival  or  exceed  fidelity 
achievable  through  simulation  and,  therefore,  represents  an  important  development 
in  model  validation.  Additionally,  generating  a  full  data  set  in  standoff  mode  for  a 
single  shot  improves  cost-effectiveness  of  funding  dollars  earmarked  for 
experimentation. 

The  fidelity  and  dimensionality  achieved  using  these  new  methods  should  prove 
useful  in  detennining  perfonnance  of  new  disruptive  energetics  and  insensitive 
munitions  and  in  determining  the  real  utility  of  fuel  additives  to  explosives.  This 
work  demonstrates  a  standoff  tool  to  characterize  explosive  performance  that  rivals 
present  modelling  capabilities,  and  should  provide  validation  data  for  physical 
chemical  explosives  models  that  incorporate  finite  rate  chemical  kinetics. 
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